Segregation-induced finger formation in granular free-surface flows

Geophysical granular flows, such as landslides, pyroclastic flows and snow avalanches, consist of particles with varying surface roughnesses or shapes that have a tendency to segregate during flow due to size differences. Such segregation leads to the formation of regions with different frictional properties, which in turn can feed back on the bulk flow. This paper introduces a well-posed depth-averaged model for these segregation-mobility feedback effects. The full segregation equation for dense granular flows is integrated through the avalanche thickness by assuming inversely graded layers with large particles above fines, and a Bagnold shear profile. The resulting large particle transport equation is then coupled to depth-averaged equations for conservation of mass and momentum, with the feedback arising through a basal friction law that is composition dependent, implying greater friction where there are more large particles. The new system of equations includes viscous terms in the momentum balance, which are derived from the -rheology for dense granular flows and represent a singular perturbation to previous models. Linear stability calculations of the steady uniform base state demonstrate the significance of these higher-order terms, which ensure that, unlike the inviscid equations, the growth rates remain bounded everywhere. The new system is therefore mathematically well posed. Two-dimensional simulations of bidisperse material propagating down an inclined plane show the development of an unstable large-rich flow front, which subsequently breaks into a series of finger-like structures, each bounded by coarse-grained lateral levees. The key properties of the fingers are independent of the grid resolution and are controlled by the physical viscosity. This process of segregation-induced finger formation is observed in laboratory experiments, and numerical computations are in qualitative agreement

A two-dimensional depth-averaged μ(I)-rheology for dense granular avalanches

Steady uniform granular chute flows are common in industry and provide an important test case for new theoretical models. This paper introduces depth-integrated viscous terms into the momentum-balance equations by extending the recent depth-averaged μ(I)-rheology for dense granular flows to two spatial dimensions, using the principle of material frame indifference or objectivity. Scaling the cross-slope coordinate on the width of the channel and the velocity on the one-dimensional steady uniform solution, we show that the steady two-dimensional downslope velocity profile is independent of scale. The only controlling parameters are the channel aspect ratio, the slope inclination angle and the frictional properties of the chute and the sidewalls. Solutions are constructed for both no-slip conditions and for a constant Coulomb friction at the walls. For narrow chutes, a pronounced parabolic-like depth-averaged downstream velocity profile develops. However, for very wide channels, the flow is almost uniform with narrow boundary layers close to the sidewalls. Both of these cases are in direct contrast to conventional inviscid avalanche models, which do not develop a cross-slope profile. Steady-state numerical solutions to the full three-dimensional μ(I)-rheology are computed using the finite element method. It is shown that these solutions are also independent of scale. For sufficiently shallow channels, the depth-averaged velocity profile computed from the full solution is in excellent agreement with the results of the depth-averaged theory. The full downstream velocity can be reconstructed from the depth-averaged theory by assuming a Bagnold-like velocity profile with depth. For wide chutes, this is very close to the results of the full three-dimensional calculation. For experimental validation, a laser profilometer and balance are used to determine the relationship between the total mass flux in the chute and the flow thickness for a range of slope angles and channel widths, and particle image velocimetry (PIV) is used to record the corresponding surface velocity profiles. The measured values are in good quantitative agreement with reconstructed solutions to the new depth-averaged theory

Les instabilités hydrodynamiques dans les écoulements granulaires géophysiques

De par le vaste champ d’applications qu’elle offre, l’étude des écoulements granulaires a connu une expansion considérable au cours de ces vingt dernières années, autant du point de vue industriel que géophysique. Ces écoulements granulaires ont une grande influence dans le monde qui nous entoure. Or, différentes instabilités hydrodynamiques peuvent naitre en leur sein, entrainant alors des changements importants des propriétés mêmes de l’écoulement. Les instabilités dues à la ségrégation par taille de particules, le développement d’ondes de surface ou encore l’apparition de ressauts dans l’écoulement en sont des exemples marquants

The kinematics of bidisperse granular roll waves

Small perturbations to a steady uniform granular chute flow can grow as the material moves downslope and develop into a series of surface waves that travel faster than the bulk flow. This roll wave instability has important implications for the mitigation of hazards due to geophysical mass flows, such as snow avalanches, debris flows and landslides, because the resulting waves tend to merge and become much deeper and more destructive than the uniform flow from which they form. Natural flows are usually highly polydisperse and their dynamics is significantly complicated by the particle size segregation that occurs within them. This study investigates the kinematics of such flows theoretically and through small-scale experiments that use a mixture of large and small glass spheres. It is shown that large particles, which segregate to the surface of the flow, are always concentrated near the crests of roll waves. There are different mechanisms for this depending on the relative speed of the waves, compared to the speed of particles at the free surface, as well as on the particle concentration. If all particles at the surface travel more slowly than the waves, the large particles become concentrated as the shock-like wavefronts pass them. This is due to a concertina-like effect in the frame of the moving wave, in which large particles move slowly backwards through the crest, but travel quickly in the troughs between the crests. If, instead, some particles on the surface travel more quickly than the wave and some move slower, then, at low concentrations, large particles can move towards the wave crest from both the forward and rearward sides. This results in isolated regions of large particles that are trapped at the crest of each wave, separated by regions where the flow is thinner and free of large particles. There is also a third regime arising when all surface particles travel faster than the waves, which has large particles present everywhere but with a sharp increase in their concentration towards the wave fronts. In all cases, the significantly enhanced large particle concentration at wave crests means that such flows in nature can be especially destructive and thus particularly hazardous

Multiple solutions for granular flow over a smooth two-dimensional bump

Geophysical granular flows, such as avalanches, debris flows, lahars and pyroclastic flows, are always strongly influenced by the basal topography that they flow over. In particular, localised bumps or obstacles can generate rapid changes in the flow thickness and velocity, or shock waves, which dissipate significant amounts of energy. Understanding how a granular material is affected by the underlying topography is therefore crucial for hazard mitigation purposes, for example to improve the design of deflecting or catching dams for snow avalanches. Moreover, the interactions with solid boundaries can also have important applications in industrial processes. In this paper, small-scale experiments are performed to investigate the flow of a granular avalanche over a two-dimensional smooth symmetrical bump. The experiments show that, depending on the initial conditions, two different steady-state regimes can be observed: either the formation of a detached jet downstream of the bump, or a shock upstream of it. The transition between the two cases can be controlled by adding varying amounts of erodible particles in front of the obstacle. A depth-averaged terrain-following avalanche theory that is formulated in curvilinear coordinates is used to model the system. The results show good agreement with the experiments for both regimes. For the case of a shock, time-dependent numerical simulations of the full system show the evolution to the equilibrium state, as well as the deposition of particles upstream of the bump when the inflow ceases. The terrain-following theory is compared to a standard depth-averaged avalanche model in an aligned Cartesian coordinate system. For this very sensitive problem, it is shown that the steady-shock regime is captured significantly better by the terrain-following avalanche model, and that the standard theory is unable to predict the take-off point of the jet. To retain the practical simplicity of using Cartesian coordinates, but have the improved predictive power of the terrain-following model, a coordinate mapping is used to transform the terrain-following equations from curvilinear to Cartesian coordinates. The terrain-following model, in Cartesian coordinates, makes identical predictions to the original curvilinear formulation, but is much simpler to implement

Solving Nonlinear Parabolic Equations by a Strongly Implicit Finite-Difference Scheme

We discuss the numerical solution of nonlinear parabolic partial differential equations, exhibiting finite speed of propagation, via a strongly implicit finite-difference scheme with formal truncation error $\mathcal{O}\left[(\Delta x)^2 + (\Delta t)^2 \right]$. Our application of interest is the spreading of viscous gravity currents in the study of which these type of differential equations arise. Viscous gravity currents are low Reynolds number (viscous forces dominate inertial forces) flow phenomena in which a dense, viscous fluid displaces a lighter (usually immiscible) fluid. The fluids may be confined by the sidewalls of a channel or propagate in an unconfined two-dimensional (or axisymmetric three-dimensional) geometry. Under the lubrication approximation, the mathematical description of the spreading of these fluids reduces to solving the so-called thin-film equation for the current's shape $h(x,t)$. To solve such nonlinear parabolic equations we propose a finite-difference scheme based on the Crank--Nicolson idea. We implement the scheme for problems involving a single spatial coordinate (i.e., two-dimensional, axisymmetric or spherically-symmetric three-dimensional currents) on an equispaced but staggered grid. We benchmark the scheme against analytical solutions and highlight its strong numerical stability by specifically considering the spreading of non-Newtonian power-law fluids in a variable-width confined channel-like geometry (a "Hele-Shaw cell") subject to a given mass conservation/balance constraint. We show that this constraint can be implemented by re-expressing it as nonlinear flux boundary conditions on the domain's endpoints. Then, we show numerically that the scheme achieves its full second-order accuracy in space and time. We also highlight through numerical simulations how the proposed scheme accurately respects the mass conservation/balance constraint.Comment: 36 pages, 9 figures, Springer book class; v2 includes improvements and corrections; to appear as a contribution in "Applied Wave Mathematics II

On dynamic interactions between body motion and fluid motion

This contribution on dynamic fluid-body interactions concentrates on applying mathematical/analytical ideas to complement direct numerical studies. The typical body may be of given shape or flexible depending on the context. In the background there are numerous real-world motivations in industry, biomedical and environmental applications, many of which involve high flow rates. A review of ideas developed over the last decade for cases of high flow rates first addresses inviscid approaches to one or more bodies free to move within a channel flow, a skimming sharp-edged body on a free surface, the sinking of a body in water and the rocking or rolling of a body on a solid surface, before moving on to more recent viscous-inviscid approaches for channel flows and boundary layers. The beginnings of certain current research projects are also outlined. These concern models of liftoff of a body from a solid surface, the impact of a smooth body during skimming and viscous-inviscid effects in the presence of more than one freely moving body. Linear and nonlinear mathematical properties as appropriate are described

A laboratory-numerical approach for modelling scale effects in dry granular slides

Granular slides are omnipresent in both natural and industrial contexts. Scale effects are changes in physical behaviour of a phenomenon at different geometric scales, such as between a laboratory experiment and a corresponding larger event observed in nature. These scale effects can be significant and can render models of small size inaccurate by underpredicting key characteristics such as ow velocity or runout distance. Although scale effects are highly relevant to granular slides due to the multiplicity of length and time scales in the flow, they are currently not well understood. A laboratory setup under Froude similarity has been developed, allowing dry granular slides to be investigated at a variety of scales, with a channel width configurable between 0.25-1.00 m. Maximum estimated grain Reynolds numbers, which quantify whether the drag force between a particle and the surrounding air act in a turbulent or viscous manner, are found in the range 102-103. A discrete element method (DEM) simulation has also been developed, validated against an axisymmetric column collapse and a granular slide experiment of Hutter and Koch (1995), before being used to model the present laboratory experiments and to examine a granular slide of significantly larger scale. This article discusses the details of this laboratory-numerical approach, with the main aim of examining scale effects related to the grain Reynolds number. Increasing dust formation with increasing scale may also exert influence on laboratory experiments. Overall, significant scale effects have been identified for characteristics such as ow velocity and runout distance in the physical experiments. While the numerical modelling shows good general agreement at the medium scale, it does not capture differences in behaviour seen at the smaller scale, highlighting the importance of physical models in capturing these scale effects

Particle-size and -density segregation in granular free-surface flows

When a mixture of particles, which differ in both their size and their density, avalanches downslope, the grains can either segregate into layers or remain mixed, dependent on the balance between particle-size and particle-density segregation. In this paper, binary mixture theory is used to generalize models for particle-size segregation to include density differences between the grains. This adds considerable complexity to the theory, since the bulk velocity is compressible and does not uncouple from the evolving concentration fields. For prescribed lateral velocities, a parabolic equation for the segregation is derived which automatically accounts for bulk compressibility. It is similar to theories for particle-size segregation, but has modified segregation and diffusion rates. For zero diffusion, the theory reduces to a quasilinear first-order hyperbolic equation that admits solutions with discontinuous shocks, expansion fans and one-sided semi-shocks. The distance for complete segregation is investigated for different inflow concentrations, particle-size segregation rates and particle-density ratios. There is a significant region of parameter space where the grains do not separate completely, but remain partially mixed at the critical concentration at which size and density segregation are in exact balance. Within this region, a particle may rise or fall dependent on the overall composition. Outside this region of parameter space, either size segregation or density segregation dominates and particles rise or fall dependent on which physical mechanism has the upper hand. Two-dimensional steady-state solutions that include particle diffusion are computed numerically using a standard Galerkin solver. These simulations show that it is possible to define a Péclet number for segregation that accounts for both size and density differences between the grains. When this Péclet number exceeds 10 the simple hyperbolic solutions provide a very useful approximation for the segregation distance and the height of rapid concentration changes in the full diffusive solution. Exact one-dimensional solutions with diffusion are derived for the steady-state far-field concentration